(will be inserted by hand later) 

Your thesaurus codes are: 

09 (09.09.1); (09.19.2); (13.25 

Detection of X-ray line emission from the shell of SNR B0540-69.3 
with XMM-Newton RGS * 

K. J. van der Heyden 1 , J. Cottam 2 , F. Paerels 2 , J. S. Kaastra 1 , and G. Branduardi-Raymont 3 

' SRON Laboratory for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, the Netherlands 

2 Columbia Astrophysics Laboratory, Columbia University, 538 West 120th Street, New York, NY 10027, USA 

3 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 

Received ; accepted 


•4); 


1 JL 1 vvyi 1V/1TX X 

AND 

ASTROPHYSICS 

October 3. 2000 ~ 


Abstract We present X-ray observations of PSR 0540-69.3 
with the XMM-Newton observatory. The spectra obtained with 
the Reflection Grating Spectrometer reveal, for the first time, 
emission from ionized species of O, Ne and Fe originating from 
the SNR shell. Analysis of the emission line spectrum allows 
us to derive estimates of the temperature, ionization timescale, 
abundances, location, and velocity of the emitting gas. 
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1. Introduction 

SNR 0540-69.3 is one of only a few Crab-like supernova rem- 
nants (SNRs) with a shell. SNR0540-69.3 harbors a pulsar, 
PSR B0540-69.3, discovered in the soft X-ray band by Seward, 
Harnden, & Helfand (1984). Like the Crab, it shows plerion 
emission from a synchrotron nebula powered by the embed- 
ded young pulsar (Chanan, Helfand & Reynolds 1984). From 
the period (50 ms) and spin down rate (4.0 x lO^^ss” 1 ) a 
a characteristic age of 1660 yr and rotational energy loss rate 
of 1.5xl0 38 ergs s“* have been derived (Seward, Hamden, & 
Helfand 1984). 

Early optical observations (Mathewson et al. 1980) classi- 
fied SNR 0540-69.3 as a young oxygen-rich SNR with an 8" 
diameter shell, bright in [O III]. Mathewson et al. (1980) also 
found a diffuse patch of [O in] ~30" to the west of the ring 
which indicated that the remnant was larger than 8". A ROSAT 
High Resolution Imager observation by Seward & Hamden 
(1994) also revealed emission from well outside the central re- 
gion, which they interpreted as from a patchy outer shell with a 
diameter of 55". This shell contributes ~20% to the measured 
flux in the ROSAT 0.1-2 keV band. More recently Gotthelf & 
Wang (2000) presented a high resolution Chandra HRC obser- 
vation of PSR B0540-69.3 which clearly shows emission from 
the outer shell. 
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In this letter, we report the results of observations of SNR 
0540-69 acquired with the Reflection Grating Spectrometer on 
board the XMM-Newton Observatory (den Herder et al. 2001). 
These observations allow us for the first time to detect and iden- 
tify emission lines from the remnant shell. Our spectral analy- 
sis allows us to derive values for the temperature, ionization 
timescale, abundances, location, and velocity of the emitting 
gas. 


2. Observations & Data Reduction 

PSR B0540 was observed repeatedly as part of the calibration 
of the Reflection Grating Spectrometer (RGS) during the Cal- 
ibration/Performance Verification (Cal/PV) phase. For the pur- 
pose of this letter, we have analyzed the 126 ks data obtained 
on 2000 May 26, which is the longest observation with the low- 
est background. The observation was performed with the tele- 
scope rolled such that the RGS dispersion axis was aligned at 
21.29°, clockwise on the sky, from celestial north. The data 
were initially processed with the XMM-Newton Science Anal- 
ysis Software (SAS). The spectra were extracted by applying 
spatial filters to the CCD image while a CCD pulse-height filter 
was applied to select the m = — 1 spectral order. The remnant 
is located on the edge of a broad region of diffuse emission 
(Wang & Helfand 1991) which complicates background deter- 
mination. However, the image is small enough on the detector 
that spatial regions offset from the source position can be used 
for background subtraction. The background spectra were ex- 
tracted by applying the same pulse height filters as used for the 
source, to spatial regions at either edge of the camera. 

A response matrix appropriate to the spatial extent of the 
source was generated as follows. A spatial mask corresponding 
to the RGS aperture was imposed on the Chandra HRC im- 
age, and the intensity distribution was integrated over the RGS 
cross-dispersion direction. The resulting profile was convolved 
with the RGS point source response matrix, generated with the 
SAS task RGSRMFGEN. The spectral analysis was performed 
using the SRON SPEX (Kaastra et al. 1996) package, which 
contains the MEKAL code (Mewe et al. 1995) for modeling 
thermal emission. 
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Fig. 1. Cross-dispersion profiles of SNR 0540-69 obtained by inte- 
grating along the dispersion axis for a: the entire spectrum and b: 
O viii Lya only. 


3. Analysis and Results 

3.J. Location of emitting region 

We initially extracted spectra with spatial filters 0.25', 0.5', 1' 
and 1.5' wide, each centered about the peak emission. The most 
interesting feature in the spectra is the presence of an Ovill 
Lya emission line that becomes more prominent as the width 
of the spatial extraction region is increased. This is an indica- 
tion that at least some of the O emission is from the SNR shell, 
and not the central region. To further constrain the spatial dis- 
tribution of the emitting gas we created a cross-dispersion im- 
age profile of the O VIII Lya line by integrating over the line 
in the dispersion direction. The cross-dispersion profile, plot- 
ted in Fig. 1, shows two peaks. Referring to the Chandra HRC 
image (Gotthelf & Wang 2000, reproduced here as Fig. 2), we 
conclude that we see O emission both from the central region 
of the remnant, centered on the pulsar, (the peak at offset zero 
in Fig. lb), as well as from the partial shell to the west (the 
second peak & 24 arcsec towards more negative offset in Fig. 
2b). 

The angular extent of this partial shell, A<j> & 30", corre- 
sponds to an effective wavelength width along the dispersion 
direction of « 0.06 A, comparable to the wavelength resolu- 
tion of the spectrometer. This implies that we do not have the 
sensitivity to resolve variations in brightness along the shell. 


3.2. Spectral analysis 

We divided the remnant into three regions. Northeast (NE), 
Central (C) and Southwest (SW), for the purpose of spec- 
tral analysis. For the NE and SW regions we use an extrac- 
tion width of 0.5'centered at ±0.375' about the central com- 
pact region, while for the brighter C region we use a width of 
0.25'. The spectra for these regions are displayed in Figs. 3 
and 4. Each spectrum shows a strong power law continuum 
with superimposed line emission, most notably Ovill Lya. 



Fig. 2. X-ray image of SNR 0540-69 .3 obtained by the Chandra HRC . 
The image contains the 50 ms pulsar PSR B0540-69.3, the plerion and 
the outer shell. The SW, C and NE regions used for spectral analysis 
are indicated. The disperion direction of the RGS is along the long 
dimension of the extraction regions. 


Emission line features are weakest in the spectrum extracted 
from the NE region, where we only detect Ovill Lya and pos- 
sibly weak OVII emission. The line emission is strongest in 
the SW spectrum where we easily identify emission lines of 
Ovill (Lya, Ly/3), Ovil (n = 2 - 1 triplet, He/3), FeXVH 
(15.01, 16.78, 17.05, 17.10 A), NeX (Lya), and Neix (n = 
2—1 triplet). Due to galactic absorption, no features can be 
seen longward of ~ 23 A. 

In the SW spectrum, O VIII Lya has a significant offset of 
A A = -0.17 A with respect to the laboratory wavelength. This 
offset must be largely due to a Doppler shift, and not to a posi- 
tional offset between the nominal pointing direction (the pulsar 
position) and the emission centroid of the shell, because the 
equivalent wavelength extent of the shell is much smaller than 
the observed shift. This wavelength shift corresponds to a ra- 
dial velocity of v ~ —2700 km s -1 . 

We synthesized a spectral model comprising three compo- 
nents: a power law to represent the synchrotron emission from 
the plerion, a non-equilibrium ionization (NEI) model for the 
thermal line emission, and foreground absorption. The free pa- 
rameters of the model are the photon index of the power law 
and its normalization, the electron temperature T e , emission 
measure, abundances, ionization age n e t of the shocked gas, 
and the column density Nu of foreground absorbing gas. Here, 
n e is the electron density of the shell, and t is the time since the 
hot gas was shocked. We adopt a distance of 5 1 kpc throughout. 

We concentrated our quantitative analysis on the SW re- 
gion since the spectrum from this region contains the strongest 
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Fig. 3. RGS 1 & RGS 2 (upper and lower plots respectively) spectrum 
of the SW region. The solid line is a model fit to the data. The gap 
between 20 to 24 A in the RGS 2 spectra is due to the failed CCD 4 on 
that instrument.The spectra exhibit a strong power law continuum due 
to synchrotron radiation from the plerion and thermal line emission 
from highly ionized atoms of O, Ne and Fe 


line emission features. Fixing the O abundance to a third solar 
(note that we have no independent way of constraining the ab- 
solute abundances, since the continuum emission is dominated 
by the plerionic emission), the abundances of Ne and Fe were 
allowed to vary with respect to O. An initial fit was made to 
the SW spectrum, using the O VIII Lya and O VI I- triplet only. 
The parameters obtained were subsequently used as starting pa- 
rameters for fits to the entire spectrum. For the Ne/O and Fe/O 
abundance ratios, we obtain Ne/O ~ 3.0 and Fe/O ~ 2.1 rel- 
ative to their solar values. The best fit model parameters are 
listed in Table 1. Fits were also made to the NE and C spectra. 
Since the line emission in these regions is extremely weak we 
fixed the NEI component parameters (i.e. electron temperature 
and ionization parameter) to the best fit parameters obtained 
from the SW region, while allowing the emission measure and 
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Fig. 4. Same as for Figure 3, but for the C region. 


photon indices of the power law component to vary. The best fit 
parameters are given in Table 1, and the spectra, together with 
the best fit models are shown in Figures 3, 4 and 5. 

As can be seen from the Figures and the \ 2 values in Table 
1, the fits are not perfect. This is probably due to a combination 
of remaining small calibration uncertainties in the RGS effec- 
tive area (up to < 20 %, shortward of 9 A), and the simplicity 
of our spectral model, which does not allow for temperature 
and ionization inhomogeneities. In particular, there is a feature 
at « 14 A, which is not fit by our isothermal model. Candidate 
emission lines are 2p - 3d in Fexvm at 14.208 A, or possibly 
2s - 3 p in FeXXI at 14.008 A, although the latter should be 
accompanied by strong emission at 12.29 A in equilibrium. In 
any case, detectable amounts of higher charge states of Fe L 
would definitely indicate the presence of hotter plasma. 


4. Discussion 

First of all, a qualitative point. We detect faint, but significant 
O vm emission from the direction of the plerion. At our spatial 
resolution, we cannot exclude the possibility that a significant 
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km s” 1 . There is evidently a wide range of ionization present 
in the plerion. Unfortunately, we do not have the sensitivity to 
try and localize the highly ionized oygen, and to constrain its 
kinematics for comparison with the [O III] image. 

Next, we check for consistency of our measured parame- 
ter values. From the apparent expansion velocity of the shell, 
v ~ 2700 km s _1 , and its size, A (j) Prf 30", we deduce an 
expansion age of 2800 yr. Comparison with the pulsar spin- 
down age then suggests that the shock has considerably slowed 
down, and that the shell has entered the adiabatic expansion 
phase. From the measured value of the ionization parameter, 
n e t & 2.5 x lO 10 cm" 3 , and the pulsar spindown age, we es- 
timate a charateristic density of n e ~ 0.5 cm" 3 in the shell. 
The value of n e t itself indicates that the shocked gas has not 
yet reached ionization equilibrium. From the measured value of 
the emission measure, we can obtain an independent estimate 
for the density. Assuming half of a spherical shell of radius 
30"= 7.4 pc, and a relative thickness of one-tenth of the radius, 
we find n e « 10 cm -3 , significantly above the value deduced 
from the ionization parameter. This discrepancy is most likely 
due to our crude estimate for the volume, and possibly by sys- 
tematic uncertainties associated with the measurement of the 
emission measure (i.e. the fact that it is based on an isothermal 
model). Alternatively, we can take it as an indication that the 
shell is definitely inhomogeneous in density, temperature, and 
ionization. 

A final consistency check can be based on calculating the 
initial supernova explosion energy from the Sedov solution for 
the shell parameters. We have, for the shell radius R s , shock 
velocity v 3 , and shock temperature T s ( e.g. Lozinskaya 1992) 

R* = 12.9<4 /5 (— ) 1/5 pc (1) 

no 


Fig. 5. Same as for Figure 3, but for the NE region. 


Table 1. Fitting results for the SW, C and NE regions of SNR 0540- 
69. 


parameter 

SW 

C 

NE 

NEI: 

n e n H V( 10 5 *cm- 3 ) 
kT c (keV) 
n c t (10 9 cm“ 3 s) 

92±6.5 

0.58±0.18 

25.22±5 

50±20 

45±18 

POW: 

photon index 

1.82±0.1 

1.2±0.04 

1.75±0.1 

ABS: 

N h ( 10 21 cm 2 ) 

3.88±1 



X 2 /d.o.f. 

2.04 

2.8 

2.9 


component of this line emission is in fact from foreground rem- 
nant gas seen in projection. If confirmed at higher spatial reso- 
lution, this detection would be interesting, because Mathewson 
et al. (1980) detected intense optical [O III] emission from this 
same inner region, with a large velocity dispersion of v « 3000 


v s = 505f 4 3 ^ 5 ( — ) 1/f5 kms 1 (2) 

no 

T„ = 3 x 10 6 <; 6/5 (— ) 2/5 K, (3) 

no 

with no = n e f 4 the density of the ambient medium, Eq the 
explosion energy in units 7.5 x 10 5 ° erg, and t± the time since 
the explosion, in units 10 4 yr. Assuming R s = 7.4 pc, v s = 
2700 km s" 1 , kT s = 0.58 keV, <4 = 0.166, and n e ~ 0.5 
cm' 3 , we find explosion energies of 2.1 x lO 50 , 1.9 x 10 51 ,and 
2.8 x 10 48 erg, derived from the estimates for the shell size, 
expansion velocity, and temperature, respectively. This range 
is too large for consistency, and the last estimate, based on the 
shock temperature, appears very low. This we take as another 
positive indication that the shell is not homogeneous. Note that 
Seward and Hamden (1994) were essentially forced to assume 
kT e ~ 5 keV in their analysis of the ROSAT HRI image, in 
order to obtain consistency with the value for the explosion 
energy estimated from the shell radius. 

5. Conclusions 

We have detected resolved X-ray thermal line emission from 
the supernova remnant shell around PSR B0540-69.3. We find 
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a spectroscopic expansion velocity of v s ~ —2700 km s _1 . 
We estimate values for the average shock temperature and ion- 
ization age, and find that the shell appears to be underionized. 
From a comparison of the geometric expansion age, and the 
pulsar spindown age we conclude that the shell has entered the 
adiabatic expansion phase. However, from apparent inconsis- 
tencies for the initial explosion energies as estimated from the 
Sedov expressions for the shell parameters, we conclude that 
the shell is inhomogeneous, which is not surprising, given its 
appearance in the high resolution Chandra HRC image. 
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